Lack of static lattice distortion in Tb2Ti207 
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Abstract. We investigated the possibility of temperature dependent lattice 
distortions in the pyrochlore compound Tb2Ti2C>7 by measuring the internal magnetic 
field distribution, using muon spin resonance, and comparing it to the susceptibility. 
The measurements are done at temperatures as low as 70 mK and external fields up to 
6 kG. We find that the evolution of the width of the field distribution can be explained 
by spin susceptibility only, thus ruling out a temperature dependent hyperfine coupling. 
We conclude that lattice deformations are absent in Tb2Ti2C>7. 
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Despite its short-range AFM correlations at temperatures lower than 100K, the 
Tb2Ti 2 07 remains in a fluctuating paramagnetic spin-liquid state down to 70mK pQ. 
This very unusual state of matter has attracted considerable attention of experimentalist 
and theorist alike. Recently it was demonstrated that, under pressure [2] high magnetic 
field [3] or both pi], Tb2Ti 2 07 does order magnetically The possibility that this order 
stems from magneto-elastic coupling was considered in all cases [21 El H]. This type 
of coupling allows the lattice to distort in order to relieve the magnetic frustration, 
concomitantly lowering the total system energy. Thus the magneto-elastic coupling, 
which is a small perturbation to the spin Hamiltonian, can select one ground state out 
of the macroscopically degenerate ground states [5j E] . In the most spectacular case the 
magneto-elastic coupling leads to long range spin order accompanied by a new lattice 
structure such as in ZnCr 2 04 [7| and CdC^CUjE]. An alternative and less dramatic 
possibility is a selection of a ground state with short range spin order and short range 
lattice deformation. In the latter case the original lattice structure is preserved on the 
average. This might be the situation in Y2M02O7 where evidence for lattice deformations 
were found by several methods [SI El [10] • 

In this work we study the possible existence of magneto-elastic coupling in Tb2Ti 2 07 
using the muon spin resonance (/iSR) technique. The basic idea is to investigate the 
nature of the changes in the local environment of the muon as the temperature decreases. 
We determine whether only the spin polarization is changing or whether the lattice 
is involved as well. Electronic spin polarization contributes to the shift of the muon 
spin rotation frequency. Lattice distortions are responsible for muon spin polarization 
relaxation. Comparing these two quantities provides information on the presence or 
absence of magneto-elastic coupling. 

Transverse [TF] and longitudinal field [LF] /xSR measurements were performed 
with powder samples on the GPS and LTF spectrometers at Paul Scherrer Institute, 
Switzerland. The measurements were carried out with the muon spin tilted by 
50° relative to the direction of the applied magnetic fields, and positron data were 
accumulated in both the forward-backward (longitudinal) and the up-down (transverse) 
directions simultaneously. This allowed us to determine both transverse and longitudinal 
muon spin relaxations rates. In Fig. [T]we show the LF [panel (a)] and the TF [panel 
(b)] data at two temperatures and applied field of 2kG. The TF data are shown in a 
reference frame rotating at a field of 1.5kG. Several aspects can be seen in the raw data: 
From the time scale it is clear that the transverse relaxation is by far greater than the 
longitudinal one. The longitudinal relaxation increases as the temperature decreases, as 
was observed previously [TT] . Finally, the transverse relaxation increases and the muon 
rotation frequency decreases upon cooling. 

The /iSR LF polarization is best described by the root exponential 



where the parameter Alp is set by taking into account the tilt of the muon spin relative 
to the longitudinal magnetic field, T\ is the longitudinal relaxation time, B^p is the 




(1) 
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Figure 1. The time dependence of the muon polarization in the (a) longitudinal and 
(b) transverse directions at an applied field of 2 kG and two different temperatures. 
The time scale is different between the two directions. The solid lines are fits to Eqs.Q] 
and [U 



background, and t is time. Similarly, the TF polarization is best fitted by a root 
exponential superimposed on a cosine oscillation 

P TF (t) = A TF exj>(-(t/T 2 )^)cos(ut + (f>)+B TF . (2) 

Here T 2 is the transverse relaxation time. The other parameters have the same mining 
as in Eq. [U The quality of the fits is presented by the solid lines in Fig. [U 

Data were collected in the temperature range 60 mK to 100 K and three fields of 2, 
4, and 6 kG. The frequency u as a function of temperature for the three different fields 
is depicted in the inset of Fig. [2j The frequency shift, 

K = (u - u)/uj (3) 

is shown in the main panel for the same temperatures and fields. We define u>o as the 
frequency of the free muon in the rotating reference frame. Similarly, we present T 2 _1 
in Fig. [3l and in the inset Tf~ . The most important aspect of the data is that all 
quantities saturate as the temperature decreases below ~2K. 

The shift and LF relaxation have a simple interpretation. K is a result of the 
sample magnetization, and T-f 1 stems from dynamic field fluctuations. The relaxation 
T 2 _1 is a bit more involved. It is a result of both static field inhomogeneities on the time 
scale of one muon spin turn, and dynamically fluctuating fields. However, since T^ 1 
is an order of magnitude lower than T 2 _1 , the contribution to the TF relaxation from 
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Figure 2. Inset: the muon rotation frequency as a function of temperature for different 
applied fields. Main figure: the shift in the rotation frequency defined in Eq. [3] 
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Figure 3. The muon spin transverse relaxation rate T 2 ~ as a function of temperature 
for three different fields. Inset, the muon spin longitudinal relaxation rates T^ 1 at the 
same temperatures and fields. 



dynamic fluctuations is negligible. Therefore, T 2 _ could be analyzed in terms of fields 
inhomogeneities only 

The T 2 _1 is related to lattice deformation via the muon coupling to its neighboring 
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spins and the system's susceptibility. To demonstrate this relation let's assume for 
simplicity that the muon is coupled only to one electronic spin S (extension to multiple 
couplings is trivial). In this case the magnetic field experienced by the muon is a sum of 
the external field H and the field from the neighboring electron H int = g/igAS where 
g is the spectroscopic splitting factor, and fig is the Bohr magneton. Here A = A(r) 
is the coupling between the muon and electron spin, which we assume depends on the 
distance between them, although at present the origin of this coupling is unknown. In 
a mean field approximation, S — > (S) = x^-Iq^b- Therefore, the muon experiences a 
magnetic field B = (1 + A%)H. Assuming that the susceptibility and the couplings are 
isotropic, the time evolution of a muon spin due to static field inhomogeneity is given 
by: 

p TF (t) = J cos[ 7m (1 + {(A)+ 5A} x )Ht]p(5A)d5A (4) 

where (A) is the mean coupling, 7^ = 85.162 MRad/ (T ■ sec) is the muon gyromagnetic 
ratio and p is the distribution of coupling variations. From Eqs. [2] and H] we find, using 
inverse Fourier transform that: (I) The frequency shift is given by 

K=^(A)X- (5) 
(II) The distribution p can be expressed as 

p(5A) = f(a/{2n\5A\})/(2\5A\) (6) 
where the effective width of the distribution, a, is given by 

a=l/\T 2X ^H\, (7) 

f(x) = yfx [sin(7rx/2) {l - 2FS(^x)) + cos(ttx/2) {l - 2FC(^x))] (8) 

and FS and FC are the Fresnel functions [12]. Thus, according to Eq. [7J if the 
ratio between T 2 l and K remains constant as the temperature decreases, the most 
likely conclusion is that a is temperature independent, and no lattice deformations are 
present. This conclusion would be wrong only if by accident (A) / (T 2 K) turns out to be 
temperature independent while {^K)' 1 does vary with T. To eliminate this possibility 
X must be measured directly as discussed below. It should be pointed out that the same 
calculation could be done by assuming that x is distributed instead of A. In this case a 
proportionality between T 2 _1 and K would mean a temperature independent distribution 
of x throughout the lattice. This would equally role out lattice distortion upon cooling. 

In Fig. H] we show (T2 / y ll H)~ l versus K with the temperature as an implicit 
parameter. The two quantities are linearly dependent at all fields. Allowing for a 
base line shift, which does not originate in localized spins will lead to a proportionality 
relation between these quantities. This stands in strong contrast to Y2M02O7 where the 
muon transverse relaxation grows as a function of x faster than exponentially. This is 
demonstrated in the inset of Fig. H] using data from Ref. [10] on a semi log scale. We 
conclude that in Tb2Ti 2 0y the muon transverse relaxation has the same temperature 
dependence as the shift. In fact, by calculating a for each data point using Eq. [7] we 
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Figure 4. Demonstrating the linear relation between the muon spin relaxation rate 
normalized by the field (T 2 7 M -ff) -1 versus the shift in muon spin rotation frequency 
K, which is proportional to the susceptibility \- The inset shows the same type of 
measurement on a semi log scale in the pyrochlore compound Y2M02O7 taken from 
Ref. [10]. 

find that Aa/a is 15% for TT^T^Or and 115% for Y2M02O7 where a and Aa are the 
average and standard deviation of a respectively. Aa/a is a measure of the relative 
change in the distances variations due to temperature changes. 

We also performed DC-susceptibility measurements using a cryogenic SQUID 
magnetometer, and home built Faraday balance [13] in a Dilution Refrigerator (DR). In 
the Faraday balance the powder sample is attached to a metallic sample holder, which 
serves as one of the plates of a capacitor, and is suspended on a spring. The sample 
experience both a field H and a field gradient of 1 T/m. The capacitance C is measured 

function of H and x °c dC 1 jdH. 

In the inset of Fig. [5] we show raw data of C _1 versus H at 70mK and 2.6K. 
At all fields the slope at 70 mK is smaller than at 2.6 K. This means that at some 
temperature, in between these two, the susceptibility starts to decrease. In Fig. [5] we 
show the temperature dependent susceptibility by both experimental methods. The 
two measurements do not overlap in the 2-5 K region. We attribute this problem to 
difficulties in cooling of the entire volume of the sample in the DR. In other words, 
the mean temperature of the sample T s is higher than the DR temperature Tdr. 
Nevertheless, as T DR decreases, and so does T s , the susceptibility does not increase. 
In fact, below Tdr ~ 0.2K, the susceptibility decreases. A saturated susceptibility and 
the decrease below ~ 0.2 K is in rough agreement with previous measurements of Luo 
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Figure 5. The DC-susceptibility as measured by SQUID, and Faraday balance 
mounted in a dilution refrigerator. The two measurements do not overlap probably due 
to the difficulty of cooling the large amount of powder in the balance. Nevertheless 
a decrease in x is observed as the temperature decreases. Inset shows raw inverse 
capacitance versus applied field. The slope is proportional to \- 



et al. |14j . This result shows qualitatively that T^" 1 and x do not vary substantially 
with T below 2 K, which supports our previous conclusions. 

To summarize, we compare the transverse relaxation rate resulting from internal 
field distribution to susceptibility measured by the shift in the muon rotation frequency 
and DC susceptibility. We find that the relaxation rate has the same temperature 
dependence as the susceptibility. This indicates that the only reason for increasing 
relaxation upon cooling is an increase in the electronic moment size. Therefore, there is 
no evidence for lattice deformation in Tb2Ti 2 07 that is static on the time scale of 0.1 
l^sec. 
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